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ABSTRACT 


This thesis is an extension of the more comprehensive 
study mede by Kenneth W. Goff (5) of the applicability of 
correlation techniques to the fleld of acoustical measure- 
ments, 


The ability of correlation techniques to separate an 
acoustical signal at a point into components according to 
(1) the source (2) the transit time from the source to the 
point in question and (3) the frequency; was used to 
measure the transient response of a reverberant room. The 
results of this correlation analysis were compared with 
pulse measurement data taken under similar experimental 
conditions. The results of the two methods were found to 
aerec. 


A technique for measuring the directive properties of 
the weve fronts associated with sound fields was devised 
and tested. It was shown that the property measured cor- 
responded with the definition of "diffusion.* Results of 
diffusion measurements taken in a reverberant room with 
various portions of the wall covered with a highly 
absorbent Fiberglas curtain were presented together with 
correspondins transient response data. 


The above mentioned experimental measurements were 
paralleled by an explanation of the theory and a sample 
of the mothematics involved in the correlation analysis 
of bounded sowmd fields. 


Thesis Supervisor: Richerd H. Bolt 
Title: Professor of Electrical Mmgincering 
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L INTRODUCTION 
1 ok HI | TOR . 


The word “correlation” implios the act of seeking 
relationships between two disercet quantities or functions 
of time or mor® properly, two onsembles of events. The 
concent is elementary. For that reason, it 16 impossible 
to pin-point some example and elite 1t as being tho first 
instance of the application of correlation. It is suf- 
ficient te say that the process is one which owes its 
development to man's imate cumlosity coupled with his 
ability to reason. 

Wistorically, correlation is a very important pro- 
ceocs. When man became dissatisfied with the nere 
recording of events and first attempted to oxnlain then 
in terms of a preceding series of incidents, he was re- 
sorting to ea form of correlation. 

Correlation, like many cther protessens involving 
basie concepts, lends itsolf to concise and explicit 
expression by methenatieal symbols, The translating of 
the process from words to mathematical symbols more or 
less parellels the evolution of the process from an art 
into 2 scloence. 

The nathematical expressim of the correlcetim pro- 
cess was a natural outgrowth of the develomont of the 
field of statistical mothods. Before 1600, no 
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mathematical conceptions of probability were recognized. 
Guilford, in his book on psychomotric nethois* cays, 
*"Ganblers had spetculsted much concerming games of chance 
wnen 1% came time to considor their losees and seins, 
particularly their losnes. They oven attempted to inter- 
ost mathematicians in their problems, thoush with cnmall 
success, ... The seventh century saw the besinning of 
serious interest in the mathematics of chance. ... 

"Bernoulli (1654-1705) publidhed the first book to 
be entirely devoted to the stibject. Demoivre (1667-1754) 
may be credited with the discovery of the normal distrlbu- 
tion curve at about 1733. From that tine on, interest was 
aroused amonz astronomers as well an mithematiclens. By 
1812 Laplece (1749-1027) hed written whet io eosidered 
the crestest single work on probability. In it ho gave 
proof of the method of least squares. 

°Tt was Gauss (1777-1054) who demonstrated the creat 
practical value of the normal curve, showing how it an- 
plied to the Gistribution of moacurenents and to errors 
mode in selentifie observations. It wac he who devised 
the fundamental methods of computation of moans, probable 
errors, and tho like. ... 

"The application of the normal curve ond elementary 
statistical methods to blologie¢sn2 and soeial date must be 


lroforence mumbors refer to the Bibliography. 
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attributed first to Quetelet (1796-1879), royal astronomer 
to the king of Delgivm. ... J41r Francis Galton (1822- 
1911) 4n working on the problems of human heredity found 
that the normal curve and its simpler applications wre 
inedequete. He invented a number of additional statistical 
tools, among them the mothod of correlation.” 

In 1808, eleven years after Galton's first pubdlica- 
tion of the emeopt developed in terms of the ideas of 
regression lines, the tern “correlation” first appeared in 
print. Pearson, Edgeworth, and Weldon further doveloped 
the mathematies of correlation to a point where it was 
generally adopted by the statisticians in the fields of 
econonices, cond aLess, biology, estronony, meteorology and 
a few othor mnysical sciences. field, in his book on 
statisties” reviews the various applications of the theory 
of correlation. As an addendum to his treatment of the 
cubjeet, Iield compiled a bibliogrephy which is excellent 
in its seope of the subject of correlation. 

Perhaps the most important advance in the field of 
statistics has been the mathematical exposition of the 
relationships betwoen theo field of conmmmicotion engincor- 
ing and statistical treatment of time series. Tho rigid 
nathematieal proofs of Norbert Wiener” not only pit the 
process of correlation on a firm methenmetical founiation, 
but also tended to accolerate furthor advancenonts in the 
field of communication onsinesring. 
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1.2 TIME SERIES AND CORRELATION 


Wiener defines time ceries as “Secuences, diserete 
or continuous, of qavantitative data assicned to snecifie 
moments in time and studied with resnect to the stetistics 
of their distribution in time. They may be simple, in 
which case they coneist of a single numerically civen 
observation at each moment of the discrete or continuous 
base sequence; or nultiple, in which case they consist of 
a number of separete auantitics tabulated according to a 
tine common to all.” 

The anplication of statictical methods to commmica- 
tion” engincering reaulres a model of an explicit function 
of time whose provertics ean be deseribed in terms of its: 

1. mean value 

2. mean-square value 

3. power spectrum 

4, correlation function 
The function of time, besides being ecanable of deserintion 
in terms of a suffielently canplicated set of probability 
densities, must also be a stationary process; that is, one 
in whieh none of the probablility dictributias which 
Geseribo the process change with time. 

It is necessary at this point to cefine the statis- 





tical cuantities associated with the measurement of random 
functions of time such aa the thermal voltage of vacuum 
tubes. 
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The mean value: If eC(V) aV is the probability 
thet V(t) lies betwoen V and V+av_ then the 


mean velue of a noise is 
& 


<vitp = fo e(v) av 


under the restriction that 


feon avei 


The mean-sduare value is 


<v* (tp se {- ev) av 


where the notation < > means the statistical average, 
as computed from the probability densities whieh deseribe 
the onsemble. 

fhe time-averaced mean value of a noise voltage 
V(t) is 


T30 


= 
V(t) > ian de fv at 
-T 


The mean-scuare time-averaged value 1s 
T 
W(t) - lam Le iG (t) at 
z7° Or -T 
A rendom function may be ascsymed to vanish outcide 


of the time interval -7T/2 <t <7/2 41n whith caso 


V(t) = face 9 (Jantt as 


= © 





where  A(f) 4s a complex function and ecauals the 
voltage spectrum of V(t) . With the aid of 
Parseval's theoren: 
-99 
{veces = frcerae = flacenr as 
“0 lg 
we may obtain an expression for the total (finite) power 


of a rondom functicn: 
V2 


Lin 4 vi (t)at = V(t) = ee aba 


Ws mals 


where W(f) » the power cpeetrum, is defined an 


lim 2 LA(zyl* 
Te 


Te auto-correlation fumetion @,,(7) of a randon 
function is 


Olt) = WEY VEST) = (v(t) v (t-7))> 


SimZlarly, the cross-correlation fimction of two 
rendon funetions is 


Ri) = THY Va TF) = Cis (t) Ve (t-7)> 


The link between the corroiation function and the 
other statistical auantities wed to desertLbe random fime-~ 
tion has its basia in Viener's theoren: the auto correls- 





aniencsity spectrum. 


The general relationshinps between the fumetion of 
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time, the spectral density, the power cpectrm, and the 
correletion fumction are shown in & schematie Aaevised by 
Kemeth ¥. Goff, Fig. 1.1. 

The results obtained by the correlation analycis of 
scouctieal problems are similar in some resnects to results 
obtoined by both steady stave methods and by pulse methods. 
Fron a rudinontary knowledge of the nature of speech and 
musical sows, ome gots the intultive feeling that, since 
the sommis normally encountered in acoustical probions 
resemble mulsed wave trains, streosing the pulse-like 
asvects of correlation methods wowld yleld significant 
results. However, such a suppecition is not sufficient 
to jJuctify the noesleet of the steady state aspects of the 
correlation process, The sbility to apply a degree of 
frequency filtering in the course of camnuting the cor- 
relation fumetion has several edvantages. Paramount among 
these several adventages is the fact thet much of tho 
present body of Imowledge on seowstical problems is based 
om tosts made with steady state, freauency filtered 
noises. ‘This fact end the ecnsequent advantages gained 
by the ability to compare results with those obtained by 
more classical methods make it desirable to compromise 
the time flitering proverties of correlation methods. 
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1.5 TOM FILTERING VS rrReQuaicyr PILPSRING 

Meany of the anplicatians of correlation teenniques 
in the f1e1d of acoustics depend unen the ability of the 
nothod to semrate signols into commonents aceording to 
their transit time from a souree. 

Fer instence, suppose two microphones were placed in 
tho sow) field of a point source radiating random noise 
{nto free space. If the microphone outputs wero mined, 
and thic combined voltage were crose-correlated with the 
Loutispeaker input voltage, the following effeetse would be 
noticed: At a delay time corresponding to the transit 
time fron the loudspeaker to the micronhone clonest to 
the source, the correlation fumetion would exhibit 2 reak 
value. At a dolay time corresponding to the trancit time 
from the loulsneaker to the second nicronhone, another 
peak would appeor in the plot of correlation fmmetion vs 
time delay. - 

If the mlerophones wore plseed too close tosether, 
the resultant plot of correlation fumetion vs time delay 
would be Gistorted. The smellest difference in trensit 
timo which produces an wiilistorted maim peak in the cor- 
relation fumction cerresvonis to the resolvins power of 
the systen. This minimum trensit time may be of the order 
ef one period of the seometrie mean frequency of the band- 
width of the random noise used; and the corresponding 
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resolvins power of the system is equal to the product of 
this minimum transit time ari the creed of sound in the 
air, 

Since the auto-correlation function and the power 
spectmm of a stationary random time function are Fourier 
transform pairs, they exhibit certain characteristic in- 
verse sproading properties. For example, a frequoncy 
spectrum, F(@) » flat from 0 to : t/r hes as ites 
transform, a function of time, f(t) , having the general 


shape of Jj sin r t vhnich has ome predominant peak and 
1 Zt 
axis crossings at $7, 2at, Tar, .., » Fig. 1.2 


(a2). As the baniwidth of F(m) AIncreases, the axis 
erossings of f(t) occur at progressively smaller values 
of time. ‘The limiting case i5 an F(m) equal tol. ‘The 
Fourler transform of such a spectrum is the so called 
"4Impulse"” or “de/fa fumetion" which is a fumetion of time 
having an infinitely small duration, an infinite ampli- 
tufe and an area of 1 umider the curve; Fig. 1.2(b). Con- 
versely, as the bandwidth decreases, the peaks in the 
function as tend to equalise., In the limiting 


ease, Fla) is composed of one freauency (= ao) , 


and the inverse FPourler trensform f(t) is 


cosine (eo) % 5; Fiss 2.2(¢). 





f(t) F(a) 


-2t -TF o BP 2r 
-E 4 


Figure 1.2(a) Fourier transform pair for finite band- 
width noise. 


f(t) F(w) 
Ox . 4 
eS 
ae) J ae 
2) ‘@) 


Figure 1.2(b) Fourier transform pair for infinitely ° 4 
wide bandwidth noise. 


F(w) 
r(t) 
oF sae 
-UL Oo A, 





Figure 1.2(c) Fourier transform pair for a pure tone. 
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From the cbove discussion end illustration it is 
obvious that the prerequisites for high resolutim in 
time are diametrically opposed to the requirements for 
high frecuency resolution. It is equnily clear that most 
acoustical meacurements made by correlatia methods ro- 
duire a compranise between the frecuency and tine 
filtering proverties of the process. 

For instance, in measuring transmission Loss of 
panels, ome soneralliy desires to mow how this property 
varies with frequency. Thus tho bandwidth of the test 
sicnal must be ewt to a degree commensurate with the fre- 
queney resolution desired. This increase in frequency 
resolution must not be so great as to create confusion 
between the successive peaks caused by the signal trans- 
mitted directly throuch the panel and those peaks caused 
by slenals arriving 2t the test voint by different paths. 
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IT INSTRUMENTATION 
2.1 IRTRODUCTI GH 


The subjoct of instrumentation must include ¢em- 
sideration of tho methods of (a) gathering data, (b) 
reducing date, (ce) presenting and recording results. 

Equipment for cathering data may include such itens 
asp electronic noise voltage Generators, loudspeakers, 
microvhones, and in a cease where the reducing equipnent 
is not portable, recording instmments. A reducing systen 
is generally composed of a number of semrate electronic 
wits each of which performs a specific process. Division 
of the reducing process into a set of sub-procosses icenis 
versatility to the system and also facilitates testing of 
the components. The end product of the reduction process 
may be presented in various ways. If the reduction system 
is electronic, the results may be presented by means of 
moter roadines ar oseLlleseepe treees. 

If a reducing system is already in existence, then the 
equimnent for satherins data and presenting results must 
be selected in such a wey that the best use will be mode 
of the features bulit into the reducing systen. 


2,2 THE ANALOG CORRELATOR 
As 4molied by the ecuetion defining the cross- 
correlation function 


DP a4 (7) = limi iG (t) fo (t-T) ats 
To 7 


oO 
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a device capeble of comoutines this fumction must perform 
the operations of (a) delay, (b) multiplication, and (c) 
intecration. If the functicns of time were olectrical 
voltages such as those produced by a micronhone in a 
sound field, the avplication of these voltages directly 
to the innut terminels of an clectrical analog conputer 
Gosigned to perform the required operation would ylold 
the correlation function directly as an output voltage 
which could be measured and recorded by properly selected 
instruments. 
The block diagram of such an analog camputer designed 

end developed by Kenneth ¥, Goff? is shown in Pig. 2.1. 
A photocranh of the correlator is show in Pig. 2.2. For 
the purposes of enalysis, the coanputer is conpoced of 
three soparate units 

1, time delay system 

2. muitiniier 

3. integrator 
Ree; of there units may be operated independently of the 
remaining components. 


2.2 (a) THE TIME DELAY 

Tne time delay mechanism employed in the computer 
being considered consists of an electro-mechanical device, 
Sis. 2.3. 498 shown in the block diagram, the tivo signal 
yoltases poins used are fed into separate channels each 

































, —?— Te ie 
3779S yV3Nit | aa, loose ! 
43 ee ig, SNONOMHONAS | | 

(gree ae | 
| | | 
37v9S 907 | | 
43040034 - te eel NAS | 
T3AST IINAVYD | yuoog NIS(L) @ | | 4 | 
| | | | 
(Sd9d 00d) | 
uy | yOlVINaOW --= | 
— = | YIddOHO ! 
2 ie 2 a | 
a = = oi 
Pica H's f sZeUP Hu (L-39) 7a f 
| (7)4 
| MYOMLIN la 
: ONILVYOSINI iy 
! MOLVHOSLNI | 
(1- eaves oo ee J 
ae x r 
(2-9) | 
sotoullede YOLOVYLBNS | | 


G3YVNOS 


YS IdILINW 


An | i 
oot y 42453 





-y¥30qv 





| 
| | 
| 


2 rd 
LV Se iy aa 


(Q3ALLINO SLNVLSNOD 


ALI IVNOILYOd0ud) 


SAIYO “LNOO 





YO ONiddjals 






JIEVINVA 





| Cute 
L- Lot 


| &, 
L-4 a1 


1, AaW130 
4) (4) a3xi4 


AV1430 4AWIL 
JIEVIYNVA 


ee lc le 


om @ 





‘IORpOTIOD Ol] JO UPaAIp yOOT¢] 





3179S ¥VINit 
HILIWWVITIIW 
JIHd Ved 









31V9S 9071 





YSISIIDWY 




















| 
) | 
| | 
| | 
43040934 INAS 
l2 Sd 00 | (QS3LLINO SLNVLSNOD 
| | 13A37 SIHAVH9 | yyoog Nis(L) @ | | rN | ALITWNOILYOd0Ud) 
| | | | 
| (Sd 00d) | 
' YsLLONd = | 4YOLVINGOW -~ | 
iene = = ==OCSCSC~S~S~<CT | Y¥3dd0HD | 
2 3 2 ane | 
a _ Fs ie 
yP(L-9) 7409) » f 3203 au (L-3) (2a, f 
| (4)4 
| MYOMLIN lm 
) ONILVYOSINI a auino> 
2 ) YOLVYOIINI | YO ONIdd31S 
(L - (a a ee — 
 (umeiinemmmmmmiinm (mim rT | | 
) i “, aW130 
(2-3)%4 JTAVINVA : 
| 
| re a 
| Jouswenteee YOLOVYLANS | | Bi v | 
| -y30Qav] | ee a | 
| , | 
| | iL oAW130 
> Pp heeiey ($+ 3) ! (3) 3 [(42-49)3 | 
| 
| avigad ant | 
4 AINdILINW 4 
| | 3nevieva | 


ass a SS i ei es 8 SS 


By 


(4) 4 








Tig. 2.2 Photograph of analog correlation computer 
syeten. 
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hevinsg its own veriable g21n omplifier, magnetic recori- 
ing head, magnetic recording track, playback head, erase 
head, and playback amplifier. The record head of channoi 
one is fixed, whereas the record head of channel two is 
movable. If tho transit times from the record to the 
playbeaek heeds for a point on the surface of the recording 
drum are equal for both channels, then no relative timo 
delay will bo introduced between the two input voltaces. 
If the trensit time between recor’ and playback heads is 
sreater on channel two than on chamel one, then 4 rela- 
tive time delay ecual to the difference in transit times 
49 introduced. For identiesl Inout voltages to oach 
channel, the output of charmmel one would be f(t) and the 
output of channel two would be f(t-T) where T is the 
difference in transit time. 

This relative time delay between channels one and 
two can be varied continuously from -15 milliseconds to 
190 milliseconds by moans of gear train and chain drive 
syatem povered by a low speed synchronous motor. 

The frequency response of the system is flat to 
within 2 2 db over a freauency renge from 100 cps to 
10 keeps. 


2.2 (b) THE MULTIPLIER 
The multiplier component of the correlator is a 


quarter=di fference-scuarins deviee which reduces the 
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process of multiplication to thet represented by the 


relation 


f. fe =} (f. + f2)? - (ff, - fa)* 


The process of vroducing a voltage provortional to 


the squere of (f, + fa) and (f, = fa) 4s per- 


formed by a squaring circuit consisting of two 6B8 
pentodes operated with ecrids in push-pull. For opera- 
tion in the region where the transconductance and grid 
voltage are lineally related, the plate current for one 
tube can be approximated by 


= 3 a 
1, = a4 + be, <j ce, 4 de, 


For the two tubes overating in parallel with a common 
plate load as showm in Fig. 2.4, the current through the 


conmon load is the sum of the plate currents 


4. = ii > in? = 2a + 2ce, ° 


is pa 


The outvut voltage 


- a mn 2 
G5 = By ~ 4 Fy = 8B + BOy 


if two squaring cirecults are operated in push-pull 
end. e, is made equal to (f. + fe) in one souaring 


and equal to (f1- fs) in the other squaring circuit, 


ete 
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; 
. 





} 
R 
& 
=] 
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Pie. Dghin plified sohemacic of the qsaringg fircuit. 
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then the output voltase will be proportionnl to the 
product f, fe : 


Variabllity in the porformance of the correlator 
fenerally results from unbalanced conditions in the 
multiplier clrevits. ‘when the multiplier is operating 
correctly, the output voltage should be zero when there 
is no input and wen there is only one input. These 
belaneed conditions correspond to the requirements that 
OxO = 0 and Oxl = O = 1x0. 

The balance for Ox0O = 0 is mede by adjusting the de 
bias of the four 6R8 tubes for ecual quiescent plate cur- 
rents whereas the balance for Oxl = 1x0 = 0 is made by 
adjusting the ac levels applied to each 6B8 tube. 

Much attention was focused upen the design recuire- 
ments necessary to reduce de drift to a degree that would 
permit continuous operation for periois up to two to 


three hours without adjustment” : 


2.2 (ec) THE INTEGRATOR 


Integration in the computer is accomplished by means 
of either a low pass filtering system or a stepping 
integrator. For acoustical measurements, a simple RC 
low pass filtering circuit is satisfactory. The corre- 
letor has several such circuits with RC times of 0.5, 
1.0, 2.0, 4, 8, and 16 seconds. ‘Selection of the proper 
RC time is another important factor governing optinum use 
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of tho ecuiment. 

The de outvout of tho RC cirevit is amplified and 
then nlaced ecross the input of a de-ac chonner which 
converts the Ae volterce inte a vronvortional 400 ens 
sipmal. Thio 400 eps signal is then amnlified again and 
placed across the output terminsia of the integrator. 

The 400 eps signel may be measured directly with a 
yacuum tube voltmeter or an oscilloscope, or it may be 
rectified and recorded by means of a de graphic amneter. 
A linear vrlot such as is obtained on a spraphic ammoter is 
extrenely useful sinee it represents the shave of the 
computed correlation fimetion. Information regarding the 
relative annilitude in decibels of suceessive peaks of the 
correlation curve is desired. Such information may be 
obtained by means of a logarithmic graphie level recorder 
which records the correlation curves directly in db. 
when used together, the linear and logarithmic plots give 
a continuous, detailed, visual record of the shave and 
magnitude ef the correlation curves as a function of tine 
delay. 

All of the information presented in the foregoing 
discussion of the analog correlation computer has been 
taken from the doctoral thesis of Kemeth Ww. Gerr”. 

Only those points whieh are considered necessary for a 
basic understanding of the computing process are pre- 


sented here. Por a more detailed treatment of the 





ho 
‘Ul 


design requirement and the performance characterictics 

of the computer, the reader is referred to either refer- 
ence 5 or to the recently completed, unpublished manual 
"Instruction Manual for Analog Correlator" by Kenneth W. 


Goff, Auguct 1954, M.I.7. Acoustics Laboratory. 


2.5 DATA GATHERING ACCESSORIES 


The accessories selected for gathering data are as 

follows: 

1. The noise source consisting of an electronic 
noise voltage generator, a power amplifier, and 
an Altee 72ER loudspeaker in a closed box baffle. 

2. The microphones used were Altec Lansing type 
21 BR-150 with 157-A bases and P5184 Altec 
Lansing power supplies. The response curves for 
these microphones ere flat to within = 2 db for 
a frequenoy range of from 100 eps to 10 keps. 
The voltage amplifiers used with these sicro- 
phones were wide band amplifiers. 

3. Sinee the reducing ecuinmnent was not portable, 
provision for recording deta was required. An 
Ampex model 350-2 twin track magnetic tane 
recorder (Figs. 2.5) was used for this purpose. 
The particular apperatus used had a frecuency 


response flat to within = 2 ab over a range of 








Pig. 2.5 Photograph of Ampex noaymwtle tone recorior 
with a tave loon. 





150 eps to 15 keps. This equipment was 
portable and yot capable of roproducing 2 
relatively distortionless replica of tho in- 


put voltage. 


2.4 ACCESSORIES FOR PRESENTING RESULTL 


Tae correlation function es computed by the anzlosg 
correlator is preeented in the form of 2 400 ens voltege 
ond a rectified de voltage. It is therefore possible to 
hneasure and record the correlation function by meems of a 
variety of electrical measuring instrumente. 

The accessories used for mensuring voltares are: 

lL. Medel 300 Ballantine vacuum tube voltmeter. 

2. Type 206 DuMont Cathade Ray Osecillescope. 


Permanent records of the output voltage were mede 
by means of: 
lL. Model AY linear seale Ssterline Angus Graphic 
Nillianneter., 


$7 
2, Bruel and Kjaer log scale Granhic Level Recorder. 
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Iii. TRANSIENT RESPONCE OF ROGHS 


3.1 INTRODUCTION 


Many studies have been made on the subject of sound 
fields in rooms. Generally there studies have been based 
on one of two distinct approaches to the problem of room 
acoustics: 

a) a study of the effect of rooms' boundaries on the 
sound field associated with steady state vnure tone 
signals. 

bd) a study of the effects of rooms' boundaries on the 
sowrxil field associated with transient signals. 


Therefore, in attempting to obtain a complete quanti- 
tative assessment of the effects of boumiaries on a sound 
field, ome might: 

a) excite the room with a non-directional steady state 
sinusoidal signal and map the amplitude and phase 
relations of the sound pressure for 

1. a11 frecuencies 
2, all possible observer positions in the room 
3. &@11 possible soureo positions in the room 

b) excite the room with a pvulse of extremely short 
duration and measure the relstive anplitude and 
direction of propogation of the associated wave 


fronts for 
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1, all time 
2. all observer positions 


3. all souree positions. 


If one has a11 of the data associated with either 
one of the above mentioned oxperiment2l approaches, he 
can mathematically deduce all of the information for the 
other mathod. For example, one can synthesize the pulse 
used in terms of a Fourier series of sinusoidals. Then 
for each freauenecy component one can pick the relative 
amplitude and phase angle of the sound pressure at any 
point in the room. By weishting these amplitude and phase 
measurements in accoriance with the Fourier analysis of 
the pulse and then addine the effects for ench freauency, 
the complete response of the sound field to a transient 
pulse can be deduced. 


~— a | 


Cn the other hand one can visualize a model of a 
room excited by a pulse. In such a model, the walls of 
the room would be elimineted and their offect on the sound 
field would be simulated by an infinite multidimensional 
array of image sourees, each of which is located in © cell 
of the same size and shape as the actual room. Each imege 
source is positioned so that the transit time of the sig- 
nal from the image source to the actual observer's 
position eausls the transit time for the signal arriving 
at the observer's position by a perticvlar sequence of 


reflections from the room's boundaries. 
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With this model, one can now specify that each non- 
directional source is in phase and is radiating o free 
traveling sinusoidal signal. By adding up the contribu- 
tion from each image sourco, one can map the amplitude 
and phase relations of a room excited by 2 steady state 
source, 

The above discussion assumes that the bounmilaries of 
the room are infinitely hard so that no phese shift or 
decrease in amplitude oceurs when a signal is reflected 
from the walls, This asswunaption simplifies the problen, 
but the resulting calculations do not conform with 
reality. It is necessary therefore to consider these 
phase shifts and attenustions occurring at the boundaries 
of the sown field, In the trensient treatment of sound 
fields one may account for the boundary effects by as- 
signing to each bounding surface an impulse response 
H(m). In studying the steady state behavior of bounded 
sound fields an absorption coefficiont and a phase shift 
are assicned to each wall. 

Any such complete analysis of a bounded sound field 
would be extremely tedious. A derire to be prectical 
would incline one to restrict the analysis to a degree 
comnensurate with the amount of detail necessary for the 
proper utilization of the room. In most rooms we are in- 
terested in producing an acoustical environment conducive 


to good audio communicctions., Therefore we may senerally 
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restrict our consideratio to the audible frequency 
range. Audio communiestions will senerally fell into 
ome of three categories, 
1. two way voeal commumications between individuels 
2, voeal commmicsations between one spealmr and an 
audience 
3, instrumental communications between musicians and 


on sudience 


Traditionally, each type of commmication mentioned 
above has associsted with it a certain physical arrangement 
of source and observer. If we accept these traditioal 
arrangements as being characteristic of the type of can- 
munication to be used in the room then we may further 
limit our analysis of a bounded sound field to those 
source and observer positions of practical interest. 

Por practical purposes we may further restrict our 
analysis by considering only those types of sounds char- 
acteristic of musical or vocel sovrees. Such simplifica- 
tions further restrict both the freauency composition of 
the testing signais and pulse duration and pulse inter- 
vais for trensient sisnel testing. 

Up to this point we have been discussing the two 
general anproaches to the study of room acoustics end the 
varLous simplifications that may be employed in this 
study. These two approaches have been quite thoroushly 
explored end the resulting body of data coupled with 
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subjective evaluations of the acoustical environment in 
rooms has led to the conclusion that variations in both 
the steady state choracteristies and the transient choar- 
acteristics of sound flelds account for differences in 
the acoustical quality of a roon, 

W.C. Sabine® recognized the importance of the 
transient behavior of sound fields and attempted to re- 
late it to a subjective evalustion of the quality of 
rooms. As a measure of transient response, Sabine 
defined a “reverberation time" the lencth of time for 
the mean square pressure of a suitably chosen distribu- 
tion of sound waves to diminish to one millionth of ite 
Roverberation time as defined by 








Sabine is a measure of tho rate of deeay of the sound 
pressures caused by a steady state sinusoidal signal 
which is suddenly discontinued. Various criterta based 
on reverberation timo ve the size of - room which is to 
bo used for a specific application, and reverberation 
time vs frequency have been devi sed! , At present, 
design curves of this tyne offer the most satisfactory 
avproach to the construction of rooms having good 
acoustics and the treatment of roans having poor acous-~ 
tical characteristics. However, it is recognized that 
atoustical design eriteria other than reverberation tine 
must be considered in the correct design cf rooms. For 
instance, one research croup? susgested that the first 
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20 db of decay of 2 sound in e room is of prinery in- 
portance in differentiating betwoen two rooms which have 
anproximately equal over-all reverberation tines. 

C.A. Mason and J. voir’, in the courre of their 
otudies of the acoustics of theaters by means of pulses 
or "tone bursts,” intimated that the time and amplitude 
distribution of reflected tono burcts might servo as an 
oid in evaluating the acoustical quality of rooms. 

In tho introduction to "Pulse Statistics ‘Analysis of 
Room Acoustics™® the authors state, *... These facts indi-~- 
cate that ... the first few reflections aro primarily 
responsible for certain inportant features of the scous- 
tleal character of rooms ... It should therefore be 
worthwhile to study this ‘short term’ transient responce 
by a method of images in which all of the wave properties 
of tho imase sourees can be considered (1.0. where the 
assunption of an inccherent source is not made). Further, 
4f an image arrey setisfying the boumlary conditions can 
be found, one should be able to treat this array statis- 
tically and thus obtain the lon; term average transient 
resvonse as well ..." 

As indicated by the above quotation, seme method of 
measurins tho transient rosponse of rooms is required for 
an evaluation of the effects of boundaries on e sound field. 
However, some estimate of the steady state behavior of 
rooms is equally desirable. One experimental method of 


moasuring the transient response of rooms has been 





sugcested “hich seems to combine some of the properties 
of both steady stete mothois and pulse methods”. This 


method employs correlation techniques. 


3,2 TM CONSLATION 4ARALYGLIS OF TRAMGIBIT ATSYQNIE£ 
Now that we have suggested some potentially useful 

mothois for making acoustical measurements, it would be 

wise to explore briefly the mathematics of this method. 
it can be shown directly that the correct applicea- 


tion of correlation measuring techniques to bounded sound 


fields will yield the impulse response of that bowmded 
sound field>. 


ceferring to Fig. 3.1, we assume that the nolse 


source 1s celivoring a signal having the complex snectrum 
F(a) to the two channels of the ecorrolator. If the net- 


works havo transfer functions H,(5) and MNa(S) respectively, 


then the outovt of the networks are: 
Falo) = Hi(jm) F, (o) (3.2-1) 


Fo(m) = He (jo) Fy fo) (3,2=2) 
fae cross power spectrum for F,(o) and Fe(w) is 


® .., (a) = a \Ha (30) F4@) Hel+je) Ft) | 


= P se (a) Ha(jm) Ha (-jo) (5,2+3) 
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where 


Po) 2 lim g |p, (o)] * (3. 2=h) 


is the power density spectrum of the noise source. 
Now the correlation is the Fourler cosine trenseforn 


of the power density spectrum 


Poy (F = [gate Ha(Jo) Ha(-Jo) 02°" ae 
a_i (5,2~5) 


Fron Fig, 3.1 we can see that frequency variation in 
the souree cen be considered as common terms to both chan- 
nels one and twe, This factor then can be considered as 


&® Gconctant, 1,.@. 


Ono (TY = Hag [ m(d0) He(+Jo) e 9°" ae 
- (3.2=6) 

If we now consider the syecial case where He = 1 
and we wish to auto-coerreiste the output and input of the 
system having a transfer fumction of Hi, (jo), the cross- 
correlation fumction is 


rd 


P 25 (tT) - D ano {m0409 e JOT aes - 


an Pg, ha (Ct) (3.2-7) 


Thus, the cross-correlation fumction is oxactly 
proportionsl to the imouwlse respmee of the system heving 
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the trensfor funetion Hi(jo). If the channel having the 
transfer function Ha(jo) 46 a roo, then the correlation 
funetian ie vreporticel to the impulse resvonse of the 
room. 

ve have telked a little about the nature of the cor- 
relation function; let us now consider its share for a 
correlation enelysis of the sound field of a room. For 
this rurpece, wo shall first define a few terms. 

Using Leplece Transformation we made the definition 


Pi. (wo) £ + ( f, (t) e “Jee a (3.2-8) 


f, (%) = {Ps (0 0 2% as 


a 


The Multiniication Theorem states: 


roan) { rato) Fela) e Jor dm = 


J esc) fo (tT) at (3.249) 


Now Gividing each side by 2T #where T = time, 
ami then taking the limit as T-- 


lim 


TT? st i Film) Felo) @ Jot doa = 


t 
lin i 
T° 2 Ln f.(t) fa(t-7) ab 2G, (7) 
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In the application of correlation teehniaues to the 
transient response of roans f,(t) 15 the voltage produced 
et tho loudcpeazer terminals by a wide load noise-voltage 
generator and f,(t) 15 tho voltage outrut of a microphoe 
placed at some observation place in the soud field. 

The voltage produced by the microphone is fe(t). It 
consists of the sum of f.(t) modified by the trensfer char- 
acteristic HN of the loudspeaker and microphonm, plus e 
swmantion of reflections es¢h of whieh onuals f, (%) 
delayed by 2 time equal to the path length of tranmission, 
1, divided by the speed of sound @ and modified by H and 
a trencsfer characteristi¢e B associated with Lts treneis- 
sion path. 

Tne frevuency domain representation cf such 2 voltage 
is 

Pe = FP, I Be meee ae do ia/8 


ame’ 8 4 «es (3,2622) 
The B's are eoaplex coefficients which eecount for the 
Tnase shifts at the reflecting surfaces, the attenuction in 
aaplitude associated with inverse square dimumition of 
rafintee sound intensity, and the attenuation in amplitwie 
eceurring et the reflectinge surface. 
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therelrore 


: - c 
Fe = Sy H @ ie Ba 


Vi: 4 


= FH = /2,/ © i * a 


Az 


If wo now make the definition: 


5 = power density spect = lim n_ik{ei © 
0 


then it follows that 


— 5¢ da 
(09) iia [28 e (3.214) 


whe re 
rio) & =| Pl) e J 
Now subatituting in the multiplication theorem: 


T} = lim e S 
Peek!) Bee {a [> “pe 


“J (?, = $e + ot) ao 


18.3% Cc J (ot Ox tls do (3.2-15)} 
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Vhen we correlate Fy and Fe using equation 3,2-15 
and 3ee-kL2 


? 42 ee fo [28s . ofa 12] 
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: 2 3 (PF (T~1_/e 
Pag (") = (ua [7 lag ~~ * fa 
= | (3,216) 


Ecustion 3.2-16 then gives tho cross-correlation 
@ qoiT) betreen the invout voltage to the loudspoaver 
end the output voltage of the microphone in terns of the 
powor density spectrum (5) of a noise voltage generator, 
the canbined transfer fymetion (4) of the loudspeaker 
and microphome, and the complex transfer fymetion cf the 
transmission math. This latter transfer fumetion 
[a,/ ¢ 9®n aeccrives the complex modification in the 


sound preasure amplitude spectrum caused by ite traveling 
the distance 1, imvolving any number of reflections. 

For the idealized case where: (a) source radiates 
noise of a flat spectrum from o, to we and zero outside 
of this spectrum; (b) the transfer responses of the loud- 
speaker, microphones, and walle are 1/03 4&t has been 
shown"° thet the eross-eorrelation curve for the loud- 
speaker ond microphone voltages wlll heve equal peaks 
for each point where the time deley correspomis to the 
transit tine from the souree to the microphone ani that 
the cross-correlation curve will damp out to a mall 
Value betwoen these points. ‘The rapidity with which it 
dauaps cut will depend upon the bandwidth (2 = o,) and 
the term (orn + m2}. 


 ———— 
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With the foregoing mathemotics in mind we may now 
vigaualice a room model that will be useful for correla- 
tion analysis. This model will be exactly the same as 
the model mentioned previously in section 3.1 in connec- 
tion with pulse analysis of rooms. In this model the walls 
are replaced by an array of image sources each of which is 
located at a distance equal to the length of the transnis- 
sion path of the reflected sount. The output of oach 
source may be modified in accoriance with the number of 
reficctions that oceur in the path of transmission from 
the scurce to the observer's position. 

In their discussion of pulse statisties theory, 

Bolt, Deak, and “vestervelt® developed the mthomatics of 
auch & model. In their model each image source occupies 

@ separate image replica of the physicel room. ‘These 
imate rooms or cells are designated by three mmbers 

(21, m, mn) each of whieh can take all intesral values from 
mimi infinity to plus infinity. The notation designating 
the imege sources is aleo an indieator of the mmber of 
reflections that the signal wmiergees in traveling from 
the source to the observer. 

We chall now make a few simple adaptations of this 
spaticl model to aid us in visualizing the correlation 
analysic of the problem. All of the sources shall be 
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Grouped im accordance with the number of reflections 
that occur in transit from the source to the observer. 
Thus the sources whose sound is reflected twice in trans- 
1% are designated eas second omer image sources. 

Tne number of image sources in each orier deponids 
upon the room coometry. For example, the model for a 
rectangular roam can be broken dawn into its image co- 
ponents as follows: 


SQunens NuspoL 

Actual i 

ist Order Inage 6 

and Order Image 18 

3rd Order Inage 38 
ete. 


4n outline of the valuet enclosed by the eetunl plus 

the six first order image sources for a reetanguler room 
ie shown in “ig. 5.2. The number of images sources of each 
order can be computed from the formule 


iil 
Number = 2 + 4 + 8 z kk (3.3=1) 


20 
where N 4a the order of the reflection and k = (N-1). 
For oxample, for the fourth order array, there will be 
2+ 16 + 8 (34241) = 66 imoge sources. 

Other factors of importance in correlation enalysis 
are the distance between imesge sources, and the distance 
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Fig. 322 


Volume enclosed by actual and first order 
Image sourees for a rectangular room. 
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from the image sources to the observer's position, 
These factors may be computed from formulas 9 and 11 
of reforence &, 
The adventages gained by reference to this spatial 
model of a room are threefold: 
a) the system gives a visual presentation of the 
problems encountered in space or time filtering. 
b) the results of the correlation analysis may be 
compared with those of other researchers who used 
pulse mothods of analysis. 
¢) the mathematics for pulse analysis developed by 
the use of this model may be applied directly to 
correlation analysis. 


34> EXPUCRIMENTAL MEASUREMENTS 


The experimontal measurements of the transient re- 
sponse of rooms were taken by cross correlating between 
the input voltage to a loudspeaker and the output of a 
microphore located in the room. A hard plaster roon, 
Room 10-3590-A, whieh had been previously analyzed by 
pulse methods” wae seleeted for the preliminary tects. 
Both the source end mierophone, as well as the test 
positions in the room, were selected to give results 
comparable to those obtained by Bolt, Doak, and 


vestervelt®. 
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The results for ome perticular test mm are shown 
in Fig. 3.3. Thowm also are the comparable results of 
& pulse analysis. In this rum, the source, a liodel 175 
James Lansing horn driving unit with a Model 1217-1290 
James Lensing throat and lens acsembly, was placed in 
the corner of the room facing the wall. ‘The microzvhone, 
em Altee 21-ER-150A was leeated 13 feet from the scureo 
alom; the short wall of the room. 

The pertinent dimensions are: 

Room dimensions: length Ix = 23 feet 
width Ly = 1364 feet 
height L2= 8.4 feet 

= 1.0 inches 
= 1.5 inches 
= 0.5 inches 
z= 1.0 inches 
= 13 feet 

= 1.0 inches 


Souree Position: 


Milkee Position: 


“#<¢4<s6 ws # 


The test signal was raniom noise generated by an 
electronic noise voltage senerator anid amplified by a 
Model 20-2 MeIntosh amplifier. A 1/3 octave band 
analysis of the loudspeaker roltage and the microphone 
voltage is shown in Fig. 3.4. The lowispeaker voltage 
and the microrhonme voltese were retorded simultaneously 
om separate tracks of a twin track tape recorder. In 
reducing the data, a section of the recording tape was 
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cut and meade into a tape loop as shown in Tig. 2.5. 

The micronhone voltage was flitered so that its spectrua 
was that passed by a 1/3 octave band filter with a geo- 
metric mean frecuency of 4000 ens. The filtered micro- 
phome voltase was fed into channel one of the correlator 
time doley and the unfiltered loudspeaker voltage was 
placed across channel two of the correlator. An RC in- 
tegrator tim of 1 secon’ and a time delay seanning rate 
of 1 millisecond per minute were used. The resulting 
linear pict of the cress correlation function together 
with the vlot for the comparable pulse test is shown in 
Fig. 3.3. In this fisure, each group of peaks in the 
eorrelotion fymection represents the arrival of the sound 
from an image source in much the sane way as the peaks 
4n the accompanying pulse snalysis. An accoapanylng 

log plot of the correlation was made. This plot is not 
presented because of the extreme physical length of the 
plotting tave. However, Fic. 5,5 contains a plot of the 
peak valuea of the correlatiom curve as 2 fimetion of C 
for Fig. 3.2. Plotted on the sem coordinates is a graph 
of the peak pulse pressure level as a function of distance 
traveled. The points for this second curve are taken from 
Fig. 4 (reference &), 
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—— INVERSE SQUARE CURVE 
EXPERIMENTAL POINTS: 
“10 4 ---- PULSE METHOD 
© -+- CORRELATION METHOD 


nm 
° 






RELATIVE LEVELS IN DB 
+ 


a. ‘57 7 10 
DISTANCE TRAVELED BY SOUND IN FEET 


Fige 3-5 Plot of the transient respomse of 
. Room 10+390A for pulse and correlstio 
analysis . 








%,.4 DISCUSSION OF EXPERIMENTAL RESULTS 


The data on pulse statistics analysis® eahows the 


spacing, of pulses arriving at the micronhone position 

when beth source and receiver are in the corner of the 

room. The correlation analysis of the transient response 

of the room should show compareble peaks in the correlation 
curve for time delays corresponding to these pulse spacings. 
The correlation curve, Fig. 3,3, shows aiscrete peaks at 
delay times of 12.7, 13.6, 21.2, 34.4, 37.2, and 47.5 milli- 
secoris. Reference to Fig, 14 of referonce 8 shows that 
these delay tiaes correspond roughly with the transit 

times for the varilous reflected pulses, 

The relative amplitudes of the peaks of the correla- 
tion give = quantitative measure of the combined effects 
of: 

1. Attenuation of signesl caused by inverse square 

diminuation of radiated sowmd intensity. 

2, Attermmation of signel caused and decrease in 
coherence of radiation from the various imece 
sources caused by 

2) The transfer fymction of ome reflecting 
surface for the case of non-degenerate 
discrete first order imge sources, 

b) The combined transfer fimetion of all the 
reflocting, surfaees for higher order nen- 


degenerate imece sources. 
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c} Decrease in coherence of the radiation from 
gBroups of degenerate imece sources. 
Thus, the plot of the correlation function va time delay 
is a measure of the short term transient response of the 
reom for a particular soures and receiver position. 

Using correlation analysis, it is alse possible to 
derive 2 value for reverberation time based on the transient 
responce of the room. The reverberation time simply cor- 
responds to the difference in delay time between the peal 
in the curve corresponding to tho arrival of sound from 
pouree 0, 0, 0 and a peak in the correlation curve which 
48 one-millionth of, or 60 db dowm from the peak of the 
direct scomd, This time may be found by extrapolating the 
rate of decay indicated by the short term transient response. 
Por the oxample discussed in the previous section where 
both source anf receiver wore in the corner, the correla- 
tion curve deeayed 20 db in the interval | = 10 mill1i- 
secomis tol= 67 milliseconds. ‘The reverberation time 
based on this rate of decay is 3.14 seconis. This value, 
then, ic a measure of the reverberation time for a 1/3 
octave bend signal having a geometric mean frequency 
ef 4000 eps. 

In line with the experimental pulse analysis of 
reference 8, additional transient response data corres- 
ponding to the date of Figs. 7 through 16 has been col- 
lected anc reduced by means of the analog correlator. 
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IV, MSASUREMENTS OF DIFFUSION 
A... INTRODUCTION 


In the literature cn room acoustics, much attention 
has been directed toward defining properties of sound 
fields which seen to have a perceptible and quantitutire 
effect wpm the coodness of the acoustical environment. 
One such property is “diffusion.* without delving into 
the technical, psycho-acoustic aspects of the problen, 
it would be advisable to diseuss some of the more intui- 
tively perceptible aspects of acoustics which are 
attencent to the property of diffusion, 

A sound fleld is defined as being campletely diffuse 
if it has wniiform energy density within the region com- 
sidered, and if the direction of propogation of the wave 
fronts at any arbltrery selected points are wholly ranéom 
in distribution. In discussing same of the trends in the 
study of soumd waves in roons!? Moree and Bolt say, *It 
now appears that reverberation time alone is not always a 
sufficient measure of auditorium excellence. It is 
desirable to hive a mean square pressure 4s nearly uniform 
as posaible over the seating area, It is also important 
to have at least a certain percentage of the sound reach 
the listener directly from the speaker, and less then 4 
certain percentage reach the listeners indirectly, after 
reflection from any single surface; ... * In the light 
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of the definition of diffuseness, the above quotation 
prescribes that for any sound field, there is an ontinun 
degree of diffusion which produces a good acoustieal en- 
virormment. One may reason that this optimum degree of 
diffusion will depend unon the use to which the sowud field 
is to bo edonted 1.,e. lecture hall, broadeast studio, 
cinema, music halls and more partieularly, upon tho 
aeoustical tradition of the population in the sound field. 
For instance, the population of Britain may be said to have 
Blightiy different preferences than the population of the 
United States. These preferences would probably be based 
upon the way the British have been used to heurings their 
music or lecturesa and also upon the peeullarities of the 
lenguage. The desree of diffusion in any sound field will 
in general be 2a fimetion of the location of the source, the 
shave of the boundarfies enclosing the sou fiold, and the 
Bhave, ebsorption eoefficient and distribution of the sur- 
faces in the sound field. 

As a sort of intuitive example of the offects of 
diffusion lot us visuealice a full symphony orchestra play- 
ine in an open field. Generally, the conductor would be 
able to heer clearly only those instruments which vere 
near to him or whose radiation was directive and beamed 
in his direction. The musicians would be ablo te hear 
proportionally less of the total effect being produced by 
the orchestra as a group, 








Under these cirewumstances, it is probable that a 
very poor performance would be tho result of this lack 
of diffusion or mixing. On the other hand, if cay a 
small string ensemble was performing in a2 very diffuse 
broadeast studio, each member of the ensemble would be 
able to hear the combined effect of tho croup, but he 
would hear relatively little of the sound thet his par- 
ticular instrument was producing. Such a condition would 
again result in an inferior performance. With the proper 
degree of diffusion, howover, it is theorized that con- 
ditions would be optimum for the production of the desired 
acoustical effect, 

In a technical report Nuwaber B 058, serial No. 1953/ 
29, dated October 19535, the British Broadcasting Company 
diseussed the property of diffusion of sow flelds. It 
was sald that imperfect diffusion may appear in at least 
two measurable forms, (1) as a chenge in pressure with a 
change in position, the frequency of the sound remaining 
econstent; (2) as a change in pressure with changes in fre- 
quency, the positions of soures and reecciver remaining 
wuniehanged. It was further stated that imperfeet diffusion 
Andicated in the steady state measurement by variation of 
pressure with frequeney or position is shown by deviation 
in the pressure, time relations from an exponential decay 
curve. These statements are perfectly reascnables hovever, 
the fumctional relationship between the degree of diffusion 
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and the trensient recponse of a sound field is not ap- 
parent, 

Correlation techniqwes seen to offer a method of 
obtaining a moasure of the diffusion of a sound field 
whith depends direetily upem the definition of diffusion, 
Moreover, this method employa almost the same oquipaent 
and arrengements as are used for the transient response 


neesurenonte. 


4,2 CORRELATION METHOD OF DIFFUSION MEASUREMENT 


Let us visunilze the following experiment. 

Ye shall place a microphone in the sound fleld crested 
by 2 lowispeaker driven by a wide bend nolse voltese 
generator. 

Now let us auto-correlate the output voltage of the 
microphone. If the noise voltage has an infinitely wide 
band flat spectrum, then it hes been shown” that the auto- 
correlation fimction of the micronnone voltage is exactly 
the combined impulse respmse of the lowlsneaker end the 
microphone. For noise voltaces of finite bandwidth, this 
auto-correlation function is the transient responce of the 
systen to 2 pulse having the corresponding, spectral com- 
position. In goneral, this auto-correlation fiumction has 
& maximum peak for zero time celay and a first axis cros- 
sins at 2 time delay depenient upon the bandwidth cf the 





noise voltage ma - @; ani the quantity mo + m1. A 
sorties of such auto-correlations is shown in Fig. 4.1. 

As explained in Chapter I, the bandwidth and the 
first axis erossing of the auto~eorrelation fumetion vary 
inversely. That is, the wider the bandwidth, the shorter 
the €itferenece in delay times between the poak cf the cor- 
relation fumetion end the first axis crossing. It is 
possible to measure the time delay correspading to this 
first exis crossing with a high degree of aceuraty. The 
time delay register associated with the analog correla~ 
tion computer pormite readings of time delay to the tenth 
of a millisecond, It 18 poselble to interpolate between 
these regicter readings by the following method, 

a} Inte each channel of the time delay, feed a 

stenderd 10 keps signal voltage. 

b) Measure the phase difference of the outputs of 
the time delay for this 10 keps signal. 

¢) Sinee the period of the 10 keps signal is pre-~ 
cisely 1/10 of a milliseeawi, it is possible to 
hand adjust the time dolay mechanism to a value 
corresponding to any fraction of a period in 
terms of degrees of phase difference between the 
time delay outputs. 

d) With the time delay mechenism thus adjusted, one 
can now remove the 10 keps sisnel, apply the noise 
voltage, and measure the voltage output correspond- 
ine to the value of the auto-correlation curve for 
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59 
that time colay. 


Now, knowing the time celay of the axis erossing, 
one can plaee two such identical microphones in 
soumd field of a souree radiating spherical sound weves. 
If the microphones are placed a distance apart equal to 
the time delay to the first axis crossins of the auto- 
correlation curve, and then orlentated in such a way that 
the soumd fronts from the source strike both mikes et the 
name instant (Fic. 4.2); then the eross-correlaticon 
fumetion at zero time delay for the output voltage of 
the two microphones will correspemd to the value of the 
auto-correlation curve for = 0. However, if the source 
microphone arrey is turned 90 degrees in 2 horizontal 
plene so that the sound waves strike one microphone, 1s 
subjected to a physical delay corresponding to the exis 
(Fis. 4.3), then the correspondins cross-correlation 
fumetion will be zero. Thus it is possible to obtain a 
plot of cross-correlation vs the orlentation the nmicro- 
phones array in degrees from 0 to 90°. If we specify 
that sero degrees corresponis to the orientation when 
poth mikes are equidistant from the souree, then we can 
convert our previously measured accurate plot of the 
auto-correlation curve vse time delay to a nist cf cross- 
correlation fumetion vs orientation in degrees by means 
of the formula =d sinc where d is the mike 
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Me, 4.62 Fhot nh of microvhone array oriented 
at O derrees in sound f40el1d of louispealmr. 
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Fig. 4.3 Photograph of microphone array orlented 
at 90 degrees in sound field of loudspeaker. 
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separation in feet, ¢ is the speed of sound in feet per 
millisecond, and is time delay in milliseconds. These 
calculations have been plotted for a noise spectrum flat fron 
100 eps to 10 keps, having a roll off of 6 db per octave 
beyond these frequencies (Fig. 4.4). 

The microphone array together with the multiplying 
and integrating cireuits of the analog correlator comprise 
an instrument which gives 2 measure of the properties in- 
volved in the definition of diffusion. 

Suppese the microphone array was placed in a perfectly 
diffuse field. Then, no matter which way the array was 
turned, the value of the cross-correlation functio would 
be the same, since by definition, a diffuse field is one 
in which the wave fronts appear to be coming wmiformly 
from all directims. It is clear that for varlations in 
Giffusion, the related cross-correlation curve would have 


some characteristic shape. 


4.3 EXPERIMENTAL RESULTS 

EXPERIMERT 2 

Two Altec 21-BR-150A microphones were pleced two 
inches apart in anechoic space in the sound field pro- 
duced by an Altee 728-B closed-box-beffled loudspeaker. 
The cross correlation of the clipped output voltage of 
the microphones was plotted as a function of orientation, 
Fig. 4.5, from 0 degrees, Fig. 4.2, to 90 degrees, 
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Fig. 4.3. The curve in Fig. 4.5 18 a measure of the dif- 
fusion of anecholic space. As should be expected, the 
source shows that the sound field produced by the loud- 
speaker in anechoic synace is not diffuse. 

The value of microphone separation used in this first 
experiment wes obtained by trial and error. It is impos- 
Bible to get the desired results by spacing the micro- 
phones the theoretically correct distance apart because 
of the finite size of the microphone diaphragms. 


EXPERIMENT 2 

The following experimental work was performed in the 
model studio, Room 20-F-O09A, located in the Acoustics 
Laboratory at M.I.T. This studio was desisned as a 
diffuse room. A physical description of this room is 
contained in reference 12. 

With the loudspeaker facing the wall showm in Figs. 
4.2 and 4.3, diffusion data was taken with the microphone 
array placed at various distances for a line perpendicular 
to the radiating face of the loudspeaker enclosure. The 
noise voltage and other equipment was the same as in 
Experiment 1. As in Experiment 1, the microphone output 
voltages were clipped 20 db below their rms value. 

The results of the diffusion measurements together 
with the plot of the transient response of the studio are 
shown in Figs. 4,6 and 4.7 respectively. 
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EXPERIMENT 3 
The following oxperimental data was taken in Room 
10-390-A, the hard plaster room. In these experiments 
the data was tape recorded and then brought into the 
Laboratory and reduced, 
The same equinment used in Experiments 1 am 2 was 
used in the following tests. 
Diffusion data and transient resnvonse deta wac taken 
fer various conditions of the room. For these tests, the 
loudspesker was plaeed facing a corner of tho room and the 
diffusion data was taken for yvarlous positions in the room. 
The positions were located at distances of 20, 15, 10, 8, 
and 5 feet fron the loudspeaker on a line between the 
loudspeaker location and the diagonally opposite corner 
of the room. The transient response data was teken at the 
15 foot nositio. | 
The above described data was taken, reduced, and 
plotted for the following room conditions: 
a) All walls of the room were bare. Figs. 4.8 and 
ALO. 

b) Cne long wall of the room was covered with a 
hishly ebsorbent FPliberglas filled quilt. Fiss. 
4.10 and 4.11. 

c) Ome long wall and one short wall were covered with 
quilting. Figs. 4.12 and 4.15. 

a) Two short walls and one long wall were covered with 
quilting. Figs. 4.14 and 4.15. 
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4,4 DISCUSSION OF EXPERIMENTAL RESULTS 

In reviewing the experimental results of Section 
4.3 1t anvears that the measuring system and the tech- 
niques employed yield results similar to those obtainable 
with a highly directional microphone. It appears thet 
the practice of setting the microphones 2a certain dis- 
tanee epart corresponds to making a micropheone which is 
extremely sensitive in measuring the direction of prop- 
ogation for ea soumd having a particular specular 
composition. Generelly, setting the microphones further 
apart makes the array more sensitive to spectra having 
lower sgoometric mean frequencies. Insufficient data is 
available to permit an evalvetion of this method of 
measuring diffusion. However, the author believes that 
the data presented herewlth is sufficient to indicate a 





possible experimental method of obtaining a measuremont 
of the degree of diffusion of a soumd field. Moreover, 
it is believed that this method warrants further study. 


4.5 SUGGESTIONS FOR ADDITIONAL WORK 

On the basis of the discussion and measurements of 
diffusion made in comection with this project, it is 
felt that there is some potential merit in the techniques 
employed. As a suggestion for additional work, it is pro- 
posed that some consideration be given to further diffusiam 
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measurements made with various bandwidths of noise voltage. 
Sueh a pnroject would make extensive experiments in some 
room such as the model studio. One possible series of 
exveriments would involve measurements for various micro- 
phone separations (selected in a manner similar to that 

of Experiment 1, Section 4.3) and various conditions of 
diffusion (adjustments in diffusion to be made by 
rearranging varlable diffusers in the model studio). 
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Ve. GCNCLUSION 


5.1 SUMMARY OF CORRELATION MEASUREMENTS 

Correlation measurements as studied in this project 
| have been of two types -~- (1) measurements of tho transient 
response of rooms and (2) measurements based on the 
definition of a diffuse sound fleld. It is an interest- 
ing coincidence that both types of measurements can be 
made with prectically the seme experimental setup. 
Transient response data may be taken by cross correlating 
between the microphone output voltages of two microphones, 
one of which is located near the source and tho other in 
the reverberant field of the room. Diffusion measurements 
are mede in the same way, by cross correlating between the 
microchone output voltages of two microphones which are 
spaced 2 particular distance apart and oriented in a 
certain way with respect to the source. 

The above illustrated coincidence serves to point out 
the versatility of correlation techniques. Unfortunsiely, 
this coincidences also serves to illustrate one of a mumber 
of pitfalls that one may encounter in the use of such 
techniques. Imvropver microphone placement in either type 
of measurement would lead to spurious results. 

It was stated before that diffusion measurements by 
correlation techniques seemed to yield results similar to 
those obtainable by use of a highly directional ateropmone. 
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It should be of interest to point out that spacings 
micronhones 2a distance apart which is calculcted from 
auto-correlation data corresponds to calibrating the 
micronvhone array for maximum directionality for the 
particular frequency spectrum of noise volitege being 
used, 

This directive microphone array constructed by cor- 
rolation techniques appears to have two major lobes in 
its directivity pattem. These major lobes are 180 
degrees apart. The calibration process consists of 
spacing the micrephones so as to eliminate any minor 
Lobes in the directivity characteristics of the array. 

It should be pointed out that the directivity which 
we haye been talking about refers only to the ability of 
the microphone array and the associated correlation equip- 
nent to discriminate between the presence of a plane wave 
front which reaches both micronhones at the same instant 
and a@ plane wave front which reaches each microphone eat a 
different instant. 

In a previously mentioned report by the British 
Broadcasting Company, 1t was concluded on the basis of 
their experimental studies that the methods of short pulse 
analysis are not satisfactory for the investigation of 
aiffusion in full seale rooms. The pulse methods refer- 
red to in this conclusion involved - (a) studies of the 


irregularity of the envelope of the decay for a short 
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pulse of sound and (b) variation in these short pulse 
decay irregularities with vosition in the room. 


15 using a direc- 


Another group of investigators 
tional microphone technique attempted to obtain a steady 
state measure of the diffuseness of sound fields. The 
conclusion to their research cited the necessity for a 
microphone that could be made highly directional over 
the whole cudio frequency renge. 

It anpvears that the correlation method of measuring 
transient response is suitable for use in full scale 
rooms. Unlike short pulse methods, correlation analysis 
of transient response may be made to yield short term 
decay rates for various pandwidths of frequencies. On 
the basis of this latter fact, it may be coneluded that 
correlation analysis ocwente a eommecting Link be tween 
room resvonse data gathered by steady state methods and 
that gathered by pulse methods. 

in addition, it ean be seen that correlation tech- 
niques, when used in conjunction with a microphone array 
calibrated for the purpose, can be made to yield data 
similar to that obtained by dlreectional microphone 
methods. Further, the microphone array can be calibrated 
for any bandwidth by use of data taken from its auto- 
correlation curve. 

From these facts it may be further concluded that the 
correlation method of analysis of bounded sound fields is 
potentially more informative than either steady state 
methods or pulse methods teken alone. 
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Pulse Statistics Analysis of Room Acoustics* 


Ro. Ti. Bott, .R. ©. DOAK, AND VP. J. WESTERVELT 
Acoustics Laboratory, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received February 2, 1950) 


Many sounds of speech and music more nearly resemble pulsed 
Wave trains than abruptly terminated continuous sounds as used 
in reverberation measurement. It is therefore not surprising to 
find that two rooms can differ markedly in acoustical quality 
even if they appear identical under reverberation analysis which 
ignores details of short transients. 

This paper introduces a pulse statistics point of view which 
takes immediate account of the pulse-like nature of common 
sounds. Fundamentally, the methed consists in examining the 
response of the room to a short pulse. The walls are replaced 
by an array of image sources (simplc images if the walls are hard, 
Or appropriately modified if there is ahsorption). These image 
arrays are then considered statistically. 

From this approach one can derive such classical quantities 
as reverberation time and mean free path. One can also analyze 


INTRODUCTION 


T is well known that reverberation time is not a 
completely adequate index of the acoustical quality 
of a room. Implicit in the concept of reverberation time 
is the assumption that the room reaches a steady-state 
diffuse condition before the source of sound is abruptly 
terminated. In practice, however, most sounds of 
speech and music can be generally classified as pulsed 
wave trains whose amplitudes and frequency compo- 
nents fluctuate sufficiently within time intervals 
shorter than the time constant of the room, so that the 
room seldom reaches steady state. Thus it would seem 
that the response of the room to transient sounds of this 
general type is an especially important physical prob- 
lem of room acoustics. The results obtained recently 
by Mason and Moir! and others? who have used short 
tone bursts to investigate acoustics of auditoriums lend 
support to this point of view. 

The task of describing mathematically the response 
of a room to an arbitrary transient, and of studying 
the roles of room geometry and distribution of absorbing 
materials in this response, is extremely complicated. 
The problem can be approached from a normal mode 
point of view,’ or one can attempt to “follow” the sound 
waves around in the room as they are reflected back 
and forth from the walls. The latter approach has been 
recently investigated by Mintzer* using Laplace trans- 
form methods. 

From the transient point of view, it is desirable to 


* This work was supported in part by the ONR, Department 
of Navy, under Contract NObs 25391, Task 7. 

1C. A. Mason and J. Moir, “Acoustics of cinema auditoria,” 
J. Elec. Eng. 88, Part III, No. 3 (September, 1941). 

* British Broadcasting Corporation, Engineering Division, 
Research Department Reports B.027 and B.035. 

$P. M. Morse and R. H. Bolt, “Sound waves in rooms,” Rev. 
Mod. Phys. 16, 117 (1944). 

‘D. Mintzer, “Transient sounds in rooms,” J. Acous. Soc. Am. 
22, 341 (1950). 


the detailed nature of discrete reflections including interference 
effects, and thus obtain an average correlation betwcen room 
geometry and the character of its pulse response. 

Idealized experiments in a hard-walled rectangular room are 
employed to illustrate the essential features of this approach. 
A point source emits an exponential damped 3600-c.p.s. wave 
train of about 2 mscc. duration. The received signals are recorded 
logarithmically on an oscillograph and the system is calibrated 
for quantitative results. Several dozen discrete reflections can be 
measured and correlated with calculation. The pulses merge into 
a more or less continuous background after a time that is calcu- 
lated and confirmed experimentally. Detailed differences arise 
according to the positions of the source and microphone in the 
room. 


use the second approach. In essence, this method con- 
sists of replacing the effect of the boundaries of the 
room by an infinite array of image sources, each image 
corresponding to one of the multiple reflections of the 
original wave emitted by the source. Finding these 
images analytically is no simple matter, and only in 
very special cases will the images be “‘mirror images”’ of 
the source. The image concept has been used exten- 
sively in earlier geometric studies®® where the source 
is considered to be incoherent. Due mainly to mathe- 
matical difficulties, little use of images has been made 
as yet In wave-acoustical investigations of rooms. 

However, it has become increasingly apparent in 
recent years that the first 20 db of decay of a sound ina 
room is of primary importance in differentiating be- 
tween two rooms which have approximately equal over- 
all reverberation times. Further, as the work of Mason 
and Moir indicates,! the time and amplitude distribu- 
tions of reflected tone bursts can be correlated with the 
acoustical quality of a room. These facts indicate that 
the images relatively close to the source, i.e., the first 
few reflections, are primarily responsible for certain 
important features of the acoustical character of rooms, 
as Brillouin has observed.® It should therefore be worth 
while to study this “short term”’ transient response by a 
method of images in which all wave properties of the 
image sources can be considered (i.e., where the assump- 
tion of an incoherent source is not made). Further, if 
an image array satisfying the boundary conditions can 
be found, one should be able to treat this array sta- 
tistically and thus obtain the long term average tran- 
sient response as well. 

This paper 1s confined for the most part to a discus- 
sion of an idealized case, a hard-walled rectangular 


§ C. F. Eyring, ““Reverberation time in ‘dead’ rooms,” J. Acous. 
Soc. Am. I, 217-241 (1930). 

6 J. Brillouin, “Sur l’acoustique des salles,” Rev. d’Acoustique 
1 (September—-November, 1932). 
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room containing a simple source which emits a short 
pulse. Statistical properties of the image array are 
investigated to illustrate the image method and to 
provide a basis for future work on more general cases. 

The limited applicability of results based on a simple 
mirror image picture has already been pointed out.’ 
However, as the discussions of both Morse and Bolt’ 
and Mintzer* indicate, a specular approximation is 
allowable in cases where the walls are not too soft, and, 
in any case, one can assume that the image of a simple 
source can be described analytically by an expansion in 
spherical harmonics, which is essentially a multipole 
expansion of the image source taken about the mirror 
image point. Thus, if the walls are not fairly hard, 
neither the simple mirror image nor the specular ap- 
proximation can be used, and other suitable representa- 
tions of the images must be found. 


PULSE STATISTICS THEORY 


1. Image Space and the Time Distribution of 
Reflected Pulses 


With these restrictions in mind, let us now set up a 
working picture for pulse analysis. We shall consider a 
simple rectangular room with perfectly reflecting walls. 
The room has dimensions L,, L,, Lz, one corner being 
at the origin of a Cartesian coordinate system. 

A sharp pulse of sound is emitted from a source in 
this room. The source is assumed to be a simple one, 
with spherically symmetric radiation. This point 
source may be located by the vector: 


r,=Xi+ Vj+Zk, (1) 


as shown in Fig. 1. A point receiver of sound pressure 


aoe bee ST ICS 


is located at a position: 
re=Ui4+ Vj+- IVk. (2) 


The vector displacement of the receiver from the 
source Is: 


Associated with this sound source is an infinite array of 
image sources each occupying, alone, an image replica 
of the physical room. This image array also is illustrated 
in Fig. 1. Each image cell is designated by three numbers 
(1, m,n), and these three numbers can take all integral 
values from minus infinity to plus infinity. The cell 
(0, 0, 0) is the actual room. 

The pulse associated with the image (1, 0, 0) reflects 
once from the wall x=L,. The pulse from (2, 0, 0) re- 
flects first from the wall x=0, then from the wall x= L,. 
Obviously the absolute value of the cell number, /, for 
cells lying along the x axis gives directly the number of 
wall reflections suffered by the pulse from the image in 
question. The pulse from (1, 1, 0) reflects once from the 
wall x=Z, and once from the wall y=Z,, so that its 
number of reflections is |/|+|m|. In fact, by this 
system of designation, the total number of wall reflec- 
tions suffered by the pulse from the image (/, m, ) is 
directly: 


Nima= |t|+|m|+|n]. (4) 
The vector position of each image source is: 
Fimn= XE VmJ Tonk. (5) 


The vector position of each image with respect to the 
recelver is: 


eS Tlmn—TIr. (6) 


To evaluate these last two equations, we note first that 






Fic. 1. An image ar- 
ray, showing image cells 
and source and receiver 
locations. The cell desig- 
nation numbers are en- 
closed in the small 
boxes. (/,=23.0 ft., 
L,=13.4 ft, 2£,=8.44 
fit.) -2 
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each image cell is a mirror reflection of the cells ad- 
jacent to it. This complicates the analysis somewhat 
in that we have two kinds of symmetry with respect to 
the basic coordinate system. lortunately the cell 
designation numbers indicate directly the type of 
symmetry: even numbers, including 0 for the actual 
room, designate cells in which the source position 
duplicates that in the original room, while odd num- 
bered cells have sources at a reflected position. There- 
fore the components for Eq. (6) are given by: 


Vn=mLy+ VY 
yn=unlh,+Z 


m= A+)DL,-X 
ym=(14+m)Ly— r l, m,n odd. 
Znx=(A+n)L,—Z 


/, m,n even, 


M=IL AN | 


(7) 


In order to generalize some illustrative calculations we 
select the longest dimension of the room as a scale 
unit and introduce the following definitions: 


fee Ly=pl, L=gl, p,q<i; 
(U—X)/L= pee, (U+ X)/L= p20, 
a Y)/L= Hye, (V+ Y)/L= yoy (8) 
(W—Z)/L=pee, OV+2Z)/L= poo. 


Thus p and gq are dimension ratios and the y’s give the 
source-to-receiver displacements for both kinds of 
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symmetry. From Eq. (6) we then get: 





L/— wre It Lm p—py-]) 
Rea 
= or + or 
C+ t-prolt) UL(n+1)p—uyol} 
{ng— pe |k even 
+ or ; or +: (9) 
C(n+ 1)q— Meo Jk odd 


Some calculations from this equation are illustrated in 
Figs. 2 and 3 which will be discussed later. 

Next let us find the average number of pulses re- 
ceived up to a specified time, ¢, after the emission of the 
pulse from the original source. These pulses come from 
all of the images within a radius | Rin,| =ct. One image 
is contained in each cell of volume V=L,L,L,. The 
number of pulses, Vj, is directly given by the volume 
of the sphere out to cf divided by the volume of one cell: 


4n|Rima|? 4act 











N > , (10) 
3V 3V 
and the number of pulses received per second 1s: 
aN, 4nreP 
= (11) 
dt V 


Oo 0 0 
0 0 0 
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Fic. 2. Graph showing N>», the number of pulses arriving at the receiver up to a time ¢ as a function of the dimen- 
sionless time parameter r. Note the extreme stepped behavior of the dashed curve for a cubical] room. 
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These equations, (10) and (11), are valid for values of 
/, m, and n, large enough so that the details regarding 
positions of source and receiver in the room can be 
neglected. These equations also are illustrated in Figs. 
2 and 3. A calculation of an actual case illustrates that 
a surprisingly large number of reflections per second are 
predicted. In a room of 10,000 cu. ft., Eq. (11) gives a 
rate of 180 pulses per second at 1/100 of a second. 
There are several conditions met in practice which 
greatly reduce this number. For one thing, the floor 
is generally quite absorptive when an audience is 
present, so that half of the spherical volume just 
assumed is essentially eliminated. Also there are 
usually a large number of degeneracies or coincidences 
which reduce the effective number of pulses, as we 
shall see. 

Figures 2 and 3 show N,, the number of pulses re- 
ceived up to a time ¢, for rooms of several dimension 
ratios and various positions of source and receiver. 
N pis plotted for convenience as a function of the dimen- 
less time parameter 


ct 
(pq)3L 


where p, g, and ZL are defined in Eq. (8). The heavy 
solid line in both figures gives the average value of V, 
as computed from Eq. (10). 

In Fig. 2 the dashed curve shows the extreme de- 
generacy of a cubical room when both source and re- 
ceiver are in a corner. The pulses arrive in large groups 
because of the symmetry of the image point lattice for a 
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cubical room and the location of source and receiver. 
The dash-dot curve for a rectangular room with p=0.5 
and g=0.33 shows similar but smaller groups of coin- 
cident pulses. Both these curves start at V,=8 since 
the direct pulse and the pulses from the seven source- 
corner images reach the receiver simultaneously at 


_t=7=0. The light solid curve is for the same rectangular 


room, source and receiver now being at point such that 
coincidences are more or less “accidental.’’ We note 
that this curve follows the average curve very closely 
for r> 2.2. For 7<2.2, the curve lies consistently below 
the average curve. This is partly because the separation 
of source and receiver are such that no pulse arrives at 
the receiver until r=1.1. In detail, this initial part of 
the curve depends in a somewhat more involved way 
on the relative positions of the real and image sources 
and the receiver. 

Figure 3 shows in more detail the initial rise of the 
curve just discussed. The dash-dot and dashed curves 
show similar initial rises for a room of p=0.8 and g=0.6. 
The same general behavior is shown by these curves, 
although both reach the average curve sooner. The 
duration of this initial departure appears to decrease as 
the room becomes more nearly cubical. 


2. Some Special Cases of Coincidences 


In Eqs. (10) and (11) it was assumed that all 
images are distinct, that is, that no two pulses arrive 
at the receiver at exactly the same time. However, if 
source and receiver are both placed in certain locations, 
it is possible for the pulses to arrive in groups, as is 


-evident from Figs. 2 and 3. It is of interest to find the 
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number of distinct pulse groups reaching the receiver 
for some of these degenerate cases. In all cases it is 
assumed that the room is rectangular, hard-walled, 
and of incommensurate dimensions. 


Case | 


Consider the source to be in a corner and the receiver 
in a position incommensurate with respect to the images 
but otherwise arbitrary. If the source is close enough 
to the corner, the images will clump together in groups 
of eight, being so close together that, effectively, each 
clump of eight acts as a single distinct image. Since 
pulses reaching the receiver will arrive in groups of 
eight, the effective number .\,’ of distinct pulses arriving 


per second is: 
14a(ct)? x(ct)? 
NV, =- =——_. (13) 
8 3V 6y 








Case 2 


Here we place both source and receiver m the center 
of the room. All the images are distinct, but because of 
the symmetry arising from the position of the source 
and receiver, a pulse from an image along a negative 
coordinate axis (the origin being taken at the center of 
the room for convenience) will arrive simultaneously 
with a pulse from an image along the positive coordi- 
nate axis. A pulse from a non-axial image in the first 
quadrant of a coordinate plane will arrive simultane- 
ously with pulses from corresponding images in each 
of the other quadrants, and a pulse from an “oblique” 
image will arrive simultaneously with pulses from cor- 
responding images in the other seven octants. Thus the 
effective distinct images are all contained in one octant 
of image space. Counting of the number of distinct 
pulses arriving can then be accomplished In exactly 
the same manner as the counting of the number of 
normal modes with frequencies less than a certain 
value. The volume occupied in image space by all 
distinct image points within a radius, cf, divided by 
the volume, V=L,L,L., occupied by each image point, 
gives for V, in this case: 


T whe 5 
NV, = —(ct?+—(ct)?-++ —(ct), (14) 
6V Soe 8V 
where: 
= Oe Oy OP 
(eh HL, Let lee 
Lr=4(L.+L,+L,). 
Similar arguments for the cases when both source and 
receiver are in a corner, at the center of a wall, or at 
the center of an edge show that all these cases (except 
Case 1, of course) can be expressed by the formula: 


1 a(ct)> @ 
~~ | +—-(n2Letnylyt nl.) (ct)? 
nanynzVL 6 = 8 


jars Mylrby+ NyqebyLl.+ Jpapoeles) wh (15) 
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where 7, ny, and 7, are to be chosen as follows: 


Ne=Ny=7.=2, source and receiver in corner; 


Hr= 7 y= 2: source and receiver at center of a 
n= 1 s-cclge, etc.; 

i source and receiver at center of a 

1=72= ! yz-wall, etc. ; 


H2=Ny=nz-=1, source and receiver at center of room. 


It is clear that for a given room the formula for V,’ will 
depend upon the position of both the source and re- 
ceiver. Further, for pulses of a finite length, for actual 
sources and receivers, and for most positions of source 
and receiver, exact coincidences will be rare, and 
“almost coincidences,” with interference effects between 
the various pulses, will be the rule rather than the 
exception. 


3. Pulse Spacing Statistics 


The close similarity between the normal frequency 
lattice and the three-dimensional image source lattice 
has already been noted. The average statistical proper- 
ties of pulses bear a close resemblance to the equivalent 
properties of the normal frequencies. Therefore, simi- 
larities in the fluctuation statistical properties common 
to both normal frequencies and pulses might be 
expected. 

The frequency spacing index y has been defined as 
the mean squared ratio of actual to average normal 
frequency spaces, for a specified frequency interval. 
This index has been evaluated for a certain class of 
rectangular rooms, and it can be calculated if normal 
frequency values are available.’ We shall see next that 
the exact analog of y for pulse spacing (defined as the 
mean squarcd ratio of actual to average time intervals 
for returning pulses) can be calculated and evaluated 
in every case for which the frequency spacing index 
index can be calculated or evaluated, provided that 
source and receiver are maintained in one corner. 

We start with the expression for the dimensionless 
normal frequency of the /, m, nth mode:? 


Himn= 3 (Pg) LP+ (m/p)?+ (n/q)? J}. 


fom m=O, 1, 2->-. 


(16) 


We obtain the time of arrival of distinct pulses from 
the image sources by using Eq. (15) for the case where 
source and receiver are both in a corner. The resulting 
equation for time of arrival can be put into the dimen- 
sionless form, 


t/4= 3(pq)*LP+ (mp)?+ (ng)? J}, 


lm,n=0,1,2---, 


(17) 


where 7 is as defined in Eq. (12). 

Notice that Eq. (17), which applies only to the case 
where source and receiver lie in the same corner, differs 
from Eq. (16) for uimn only through the appearance of 
(p)—! and (q)~ in place of p and gq. Since p and gq enter 


7R. H. Bolt, ‘‘Norma! frequency spacing statistics,” J. Acous. 
Soc. Am. 19, 79 (1947). 
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symmetrically in both Eqs. (16) and (17), they can be 
interchanged. Thus in order to convert any equation 
giving the frequency spacing index into a form which 
will yield the pulse spacing index, one need only substi- 
tute p— for g, g~! for p, and 7/4 for w. Of course, if the 
room proportions are such that the required weighting 
factors cannot be evaluated, but actual pulse spacing 
values are available, for example from experiment, y 
for pulse spacing is easily evaluated from its defining 


equation: 
1 B Otinn a 
E (=) rs 


TB—TA A (OT) ay 


where (67)s,/4 is obtained by making the substitutions 
p—(q)", gq (), uw 7/4 In the equation giving the 
average spacing between adjacent normal frequencies.’ 








(18) 


Vas= 


4. A Derivation of the Mean Free Path 


Pulse analysis yields a straightforward derivation of 
the classical mean free path. We assume that the order 
numbers are so large that the coordinates of the source 
and receiver can be neglected (that is, source and re- 
ceiver can be considered to be at the origin). Using 
polar coordinates (7, 0, ) we can express the number of 
reflections associated with each cell as: 


cos@ cos@ sind cos¢ sind 
Nine] —— t+ | (19) 
Lc Ly i 


We next obtain an average value of the number of 
reflections out to a given radius r=cl, by averaging 
Eq. (19) over all angles: 


Nimrdo 
rs 


g 


Nimn= (20) 


So AV 
where do=r* cos¢d¢d6, So= (427?) /8, S=2(L,L.+- LL, 
+L,£,), and the integration is taken over one octant 
of space. By definition, the mean free path is equal to 
the total distance traveled by an average pulse in a 
given time ¢, divided by the number of reflections of 
that pulse during the same time. Therefore: 


miei ye — 41/5. AM 


5. Energy in the Pulse 


We now consider the energy contained in a sound 
pulse and follow the course of that energy as the sound 
becomes dispersed throughout the room and absorbed 
at its boundaries. The total power radiated by a simple 


source 1s :° 
us pv’Qo" 
= —((po")ro”) ergs/sec., 


II = 
2c pc 


8 P.M. Morse, Vibration and Sound (McGraw-Hill Book Com- 
pany, Inc., New York, 1948), second edition. 


(22) 
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where Qo is the source strength, and (92) is the mean 
square sound pressure at a distance rp from the source 
in a free field. This equation applies to steady-state 
radiation. It is also valid for a single pulse wave train 
which has a sufficiently narrow spectral distribution. 
This requirement is fairly well satisfied, for example, 
if there are ten or more waves in the train, and if the 
wave amplitude is fairly constant throughout the 
duration of the pulse. We designate the length of this 
pulse 7, and write the total energy contained in the 
pulse as: 


E,= Urp= (Ar/pc)((po")ro”) Tp ergs. (23) 
The energy density in the pulse is continually diminish- 
ing as the pulse radiates outward. At any instant of 
time, ¢, after emission, the volume of space containing 
the pulse is: 

V p= (42/3) ALP —(t— 1p)? |J~Arer,. (24) 
We next consider the multiplicity of pulses arriving at 
the receiver from all directions as time progresses (still 
assuming that the room is lossless). We take the number 
of pulses arriving per second, as given by Eq. (11), and 
multiply this by the energy density during the duration 
of a pulse passage, obtaining the energy density per 
second (provided the pulses add incoherently): 


los aN » 4r (po )ro? 


aa is (25) 
V, dp V 


W, 





ergs/Cc-sec. 





If we now multiply W, by the fraction of a second 
occupied by the individual pulse, we obtain the average 
energy density in a room of volume V: 


An (por)rer 
Ve=— 
pc 





T, ergs/ Ce (26) 


A more convenient quantity experimentally is the mean 
Square pressure in the room: 


4ac((po?)r0°) Tp 


(PRr)= p?PW r= (27) 


Suppose we take a pulse of +,=20 msec. duration 
which generates a mean square pressure of one dyne 
(74 db pressure level) at a distance of one meter from 
the source in a free field. If this pulse is emitted in a 
room of 10,000 cu. ft. volume, we find that the average 
mean square pressure throughout the room after the 
pulse is dispersed, as calculated from Eq. (27), is about 
69 db, or only 5 db less than the sound pressure level 
in the original pulse as measured at one meter. In 
practice, the sound is rapidly dissipated by absorption 
and may be canceled out by interference effects. 
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Fic. 4. Graphs of peak 
pulse pressure level as func- 
tions of distance from the 
source for several sequences 
of cirect and reflected 
pulses. 
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6. A Derivation of the Eyring 
Reverberation Equation 


Let us assume that all of the boundaries of the room 
are equally absorptive, and that their absorption is less 
than 20 percent, so that the specular reflection ap- 
proximation is valid. Further, let us suppose that the 
energy in a pulse is diminished upon reflection by a 
factor (1— a), where a@ is the average over all angles of 
the free wave absorption coefficient. Then the pulse 
energy associated with the .Vth image is: 


a — E,(1 = a)V= el — q) (tS /4¥) 


which leads at once to the Eyring reverberation equa- 
tion: 


(28) 


To= ETS Inf 1/(1—a@) ]. 


This is the average long term reverberation time of the 
room. To study the short term response for a fairly 
hard-walled room, one could utilize the specular reflec- 
tion approximation*‘ for each of the first few reflections 
and thus obtain in detail the approximate short term 
response of the room to a particular transient. Further 
refinements of the long term response might be obtained 
by treating groups of images statistically, just as 
modified decay equations have been obtained by group- 
ing of normal modes.* 


PULSE MEASUREMENTS IN A HARD-WALLED 
RECTANGULAR ROOM 


1. Experimental Procedure 


In order to obtain some simple results which could be 
interpreted from a pulse statistics point of view, short 
sound pulses were produced experimentally in a hard- 
plaster-walled rectangular room, and photographs of the 


sound arriving at a microphone were made using a 
cathode-ray oscillograph. The source used was a W.E. 
713A receiver unit feeding into a {-in. diameter brass 
tube 8 in. long packed with steel wool to present a high 
acoustic impedance to the diaphragm. The effective 
source was the end of the tube, which was small and 
could be considered approximately a simple source. 
The microphone was a W.E. 633A dynamic type, the 
output of which was passed through an ERPI RA-363 
octave filter and through the amplifier section of a 
ERPI RA 277-F sound analyzer. 

The output of the analyzer was attenuated logarith- 
mically by a Kay Labs Type 510-A Logaten and then 
put across the vertical detlection plates of a DuMont 
247 cathode-ray oscillograph. A pulsed carrier of about 
2 msec. duration was produced by mechanical switching 
of the speaker input which consisted of a 3600-c.p.s. 
signal from a Hewlett Packard 200 D audio oscillator 
amplified through a Fairchild audio amplifier. A slow 
(approximately 42 cm/sec. on the screen) external single 
sweep on the oscillograph was activated mechanically 
shortly before the beginning of the pulse. Thus the 
logarithm of the acoustic signal, picked up by the micro- 
phone during the 0.25-sec. interval after the pulse was 
emitted, was recorded linearly on the oscillograph 
screen as a linear function of time, and was photo- 
graphed. 


2. Calibration and Auxiliary Data 


The horizontal sweep speed was calibrated by a 
30-c.p.s. signal direct to the vertical plates from the 
oscillator. The vertical deflection was calibrated in 5-db 
steps, using the attenuator pad on the ERPI analyzer. 
The horizontal sweep was found to be linear within 
experimental error, while the vertical deflection was 


SS) 


found to be lmear in db with a maximum deviation of 
+1 db over the main working range of 50 db. With 
somewhat larger variations from linearity at low ampli- 
tudes, the vertical deflection is approximately linear 
over a range of 70 db, which includes background level 
in the pictures. 

The following assumptions and conventions were 
adopted in interpretation of the data: (a) All pictures 
were calibrated in db vs. ft., since time in the room 
corresponds to distance in image space; (b) it was as- 
sumed that the acoustic pulse emitted by the source in 
all pictures was the same in both shape and height; (c) 
a decibel reference level was established in terms of the 
peak amplitude and was taken as zero db at 0.1 it. 
from the source, all other levels thus coming out nega- 
tive ; (d) the velocity of sound was taken as 1120 ft./sec. 

In order to get a rough check on the validity of the 
assumption that the actual source was a simple source, 
two sets of pictures were taken in which source and 
receiver were placed near the center of the room (so that 
direct pulses would arrive without interference from 
images) and were spaced several different distances 
apart from 0.50 to 4.0 ft. The peak levels of these direct 
pulses plotted against distance between source and 
receiver on semi-log paper are shown in Fig. 4. The best 
straight line fit to the experimental points is a line with 
inverse square slope. This in turn indicates that in- 
terpreting the actual source as a simple point source is a 
fair approximation. Another check on this approxima- 
tion was made by measuring peak pulse amplitude in 
db at a constant radius from the source, but at six 
different angles with respect to the tube, from 0° to 
180°. These measurements showed a maximum varia- 
tion of -+2 db in the peak values. 


3. Pulse Shape and Energy 


Pulse length and carrier frequency were selected 
empirically. The pulse length was chosen short enough 
compared to room dimensions to give at least several 
clearly separated ‘‘echoes,” but long enough to include 
at least five or ten cycles of the carrier. In turn the 
carrier frequency was selected as high as possible 
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Fic. 5. Photograph for a non-degenerate case in which the 
distance between source and receiver is 1’. This photograph shows 
a typical direct pulse without interference from image pulses. 
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Fic. 6. Graphical enlargement of the direct pulse shown in Fig. 5, 
illustrating the general shape of the envelope. 


(without encountering large air absorption) in order 
that the speaker would operate in a region of high 
efficiency and good transient response. The sweep speed 
and carrier frequency were such that carrier details 
are barely unresolvable, although individual oscillations 
of the carrier can be distinguished in some places in the 
pictures. 

The exact nature and shape of the pulse was subse- 
quently determined by analysis of the pictures (e.g., 
Fig. 5). The average pulse dimensions as determined 
from the photographs are illustrated in Fig. 6, which isa 
graphical enlargement of the direct pulse shown in 
Fig. 5. The detailed shape of the pulse is represented by 
a stepped curve in which the width of most of the steps 
is equivalent to several cycles of the carrier. At the 
onset of the pulse the first two or three swings of the 
alternating and exponentially increasing carrier are just 
discernible and the pulse builds up to maximum ampli- 
tude in about two cycles of the 3600-c.p.s. carrier. The 
pulse is almost symmetric but its peak amplitude in the 
positive direction is measurably greater than in the 
negative direction for this particular picture. (Sym- 
metry of the pulse, of course, depends upon the phase 
of the carrier at the instants of switching it on and off.) 
Following the peak, the pulse decays in a rigorously 
exponential manner within experimental error as indi- 
cated in Fig. 6 by the two straight lines drawn through 


B-O-L. T , DO rai 
the stepped curves. The pulse «decay rate is about 
2.5X 10° cdb/sec., and depends principally upon the 
loudspeaker characteristics. The background noise 
showing on Tig. 5 is mainly a.c. power ripple. Maximum 
background level for all pictures was —65 db. 

The energy in a pulsed carrier having an envelope 
that rises instantaneously and then decays exponen- 
tially is 
w+ 2" 


y(w?+ y’) 


where fo is the peak pressure in dynes at rp cm, w 1s the 
angular frequency, and y is the exponential decay con- 
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stant. For the pulse shown here y=580 and w= 72007 


so that (y/w)*<1. Hence Iq. (29) reduces to 


T pyro 


(30) 
pcy 


If we let p, be the r.m.s. pressure in the room a rela- 
tively long time after emission of the pulse (i.e., when 
the sound is fairly umformly distributed), then the 


ratio p./ po is, using Eq. (27): 
(=) 
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Fics. 7-12. Series of pulse photographs showi ing response in a hard-walled rectangular room when the source is placed in 
a corner and the receiver is located at various points around the room. Note the extreme number of coincidences when the 


receiver is also in a corner (e.g., Figs. 8 and 11). 
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Fic. 13. The total 
number of distinct pulses 
N, arriving at the re- 
ceiver, plotted as a 
function of time, when 
the source and receiver 
are both in the corner 40 
of the room, asin Fig. 11. 


or, using Eq. (30) and letting 7,—3y: 


Ps (=) ; 
Po V 


This expression neglects dissipation of energy in the 
room which is easily taken into account if the reverbera- 
tion time of the room is known. 


(31) 


4. The Period of Resolution 


A period of resolution ¢, can be defined as the time at 
which the expected interval between successive pulses is 
just equal to the effective pulse duration 7». At a time 
t, after the first pulse, one would no longer expect to 
‘see’ the individual echoes in the clear, but rather a 
smear characterized by an envelope above which the 
peaks of the pulses would occasionally appear. 

An expression for the period of resolution ¢, in terms 
of the pulse width 7, can be obtained from the correct 
(i.e., appropriate for the existing number of degenera- 
cies) expression for V, or V,’ given in Eqs. (10), (13)- 
(15) by solving the difference equation for 1,: 


Neo (32) 


This has been done for the conditions of the experi- 
ments performed and leads to the following values for 
the period of resolution: 





Source and receiver in !-=9.98X 1073 sec. 


arbitrary positions: 


Source in corner, receiver tr=26.3X 10-3 sec. 


in arbitrary position : 


Source and receiver in t-=56.8X 1073 sec. 


corners: 


These values correspond to an effective pulse duration of 
Tp=1.74X10-3 sec., which is the pulse width at an 
amplitude equal to 1/e times the peak amplitude. The 
resolving time as defined is depicted in each of the 
pulse pictures by a vertical line appearing, in every case 
but one, to the right of the initial pulse. The individual 
pictures will be discussed in detail later. 


5. Envelope and Decay of Unresolved Sound 


The height of the vertical line that mdicates the 
resolving time has been adjusted to equal the r.m.s. 
sound level relative to the peak pulse level as obtained 
from Eq. (31). The subsequent decay of this more or 
less continuous sound is portrayed by the two slightly 
slopmg horizontal lines that form the envelope of the 
unresolved sound. The slopes of the envelope were de- 
termined from the measured reverberation time of the 
room at 3600 c.p.s. 

It will be noted that the size of the envelope is, in a 
number of cases, considerably too large (e.g., Figs. 8, 
10, 11). In every case where the theoretical envelope 
does not fit the photograph, the source, or both the 
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Fic. 14. The aver- 
age and actual rates 
of arrival of pulses at 
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both source and re- 
celver are in the 
corner of the room 
as in Fig. 11. 
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Fics. 15 AND 16. Pulse responses for the receiver in other locations. 


source and receiver, was in a corner position. Since the 
source was of finite dimensions it was never exactly in a 
corner but usually from 0.2 to 0.5 of a wave-length from 
one or more of the walls. Hence there was interference 
between the source and its images, resulting in a strongly 
directional effective source, and, consequently, in direc- 
tional image clusters. Thus for the source in a corner 
and a particular location of the receiver, certain image 
sources will contribute very little effective energy to 
the receiver. This amounts to saying that the source 
effectively puts less energy into the room. This point 
will be made more evident from the discussion of the 
photographs. 


6. Discussion of Pulse Photographs 


The photographs were taken in a room of dimensions 
L,=23.0, Ly=13.4 and L,=8.44 ft. To specify the 
location of the source and receiver, we shall use (X, V, Z) 
for the source position and (U, V, WV’) for the receiver 
position, distances being measured in feet. 

In Fig. 5, the source coordinates are (11, 7.2, 3.9 ft.) 
and the receiver was located at (12, 7.2, 3.9 ft.). This 
position is non-degenerate and so the resolving time 


shown is /-=9.98 msec. The envelope is seen to be 
shghtly too large. This can perhaps be accounted for 
by the fact that both receiver and source are near the 
center of the room so that images from opposite pairs 
of walls interfere with each other consistently. 

For Fig. 7, (XYZ) are (1, 1.5, 0.5 in.) and (UVW’) are 
(5.75, 3.35, 2.11 ft.). The source is thus effectively 
eightfold degenerate, being in the corner, so the resolv- 
ing time is /;= 26.3 msec. It will be noted that the direct 
pulse is exceeded in height by two succeeding pulses. 
The peak amplitudes of the pulses indicated by the 
squares are plotted on Fig. +. This sequence corresponds 
to pulses from the successive images at (0,0, 0), 
(0, 0, 2Z.), (0, 0, 42.), and (0, 0, 6L,), each of which is 
a cluster of eight image sources. The differences in 
amplitude show clearly that each image cluster has 
distinct directionality, with especially great variations 
in the XY, Y plane. The long-dashed line on Fig. 4 is the 
peak level one would expect for incoherent addition of 
16 sources (i.e., 12 db above peak for a single source). 
It is seen that the direct sound is 12 db below this line. 
However, the sound from the image sources lies along 
a line 14 db above the peak for a single source, which is, 
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in turn, 4 db below the level expected if the eight sources 
in each cluster all added coherently. 

Figures 8-12 are a series in which (YYZ) are still 
(1, 1.5, 0.5 in.) and in which U=1 in., H’=1 in. for 
all photographs. In Fig. 11, V=£,—1; in Fig. 12, 
V=—3L,/4; in Pig. 13, V=L,/2; and in Figs. 14 and 15, 
V=1 in. The directionality of the image clusters 1s 
clearly evident in Fig. 8. Here there is no distinct ‘“di- 
rect” pulse since the pulse from the image at (0, 2L,, 0) 
arrives simultaneously with the pulse from the actual 
source. Evidently these two pulses interfere since the 
height of the first pulse is small. The series of tall 
pulses is again from the images at (0, 0, 2L.), (0, 0, 4.) 
etc., each of which arrives simultaneously with corre- 
sponding pulses from the images at (0, 2L,, 2L.), 
(0, 2L,, 4£,), and their reflections in the X, Y plane. 
Each of these pulses then represents the simultaneous 
contributions of 32 image sources. As is evident by 
comparison of the peak heights with those in Fig. 7, 
the sources do not add coherently, however. The 
marked effect of interference is thus responsible for the 
low value of the actual envelope in this picture. The 
symmetry of the images corresponds closely to the sym- 
metry when both source and microphone are in the 
same corner so the resolving time indicated is ¢,=56.8 
msec. 

Figures 9 and 10 show effects similar to those Just 
discussed for the other positions of the microphone 
mentioned previously. The resolving time for these is 
again ¢,= 26.3 msec. since the microphone position is no 
longer ‘‘degenerate.” The influence of image-source 
directionality continues to be apparent. 

In Figs. 11 and 12 both source and receiver are in the 
corner. Figure 12 serves to mdicate the reproducibility 
of Fig. 11 which was possible with the experimental 
arrangements employed. The resolving time here 1s 
t,= 56.8 msec. The triangles refer to the points shown on 
Fig. 4. These pulses each represent the simultaneous 
contributions of the images at (0, 0, +2Z2,), (0,0, +42.) 
etc. As is evident from Fig. 4, the peaks very nearly 
follow the inverse square law. The peak amplitude, 
however, corresponds to 9 db above peak for a single 
source, whereas each pulse comes from 16 sources. 
The resulting peak in this case is less than would be 
expected for incoherent addition. Consequently it is 
again not surprising that the actual envelope is below 
the predicted value. The lines drawn on the lower half 
of the photograph indicate the calculated times of ar- 
rival of distinct pulses. Figure 13 shows the total num- 
ber of distinct pulses .V,’ arriving up to a time ¢/, as 
calculated from Eq. (15) for this location of source and 
receiver; and Fig. 14 shows the average rate of arrival, 
and a graphical portrayal of the actual pulses arriving. 
Each line on Fig. 11 corresponds to a pulse on Fig. 14. 
It is seen that agreement is quite good. 

In Fig. 15 the source is still in the corner but the 
receiver is at (4L:, +L,, 1 in.). This picture should be 
compared with Fig. 10. 
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In Fig. 16 the source remains at (1, 1.5, 0.5 in.) while 
the receiver is placed in the opposite corner (nearly) at 
(L,, Ly, 0). Here the pulse from the source arrives al- 
most simultaneously with the pulses from the three 
image clusters at (0, 2, 0) (2L.,, 2L,, 0), and (212, 0,0). 
The resulting interference is clearly evident from the 
appearance of the first pulse. In this case the symmetry 
is again almost the same as that for both source and 
receiver in the same comer. Therefore t-=56.8 msec. 
is used. Agreement of actual and calculated envelopes 
is seen to be fairly good. 

Although the present experiments are restricted to an 
idealized room, a close connection can already be seen 
between these results in a small room and pulse studies 
in large auditoriums. For example, in Fig. 17 the source 
is at (1.5, 6.7, 1 ft.) and the receiver at (20,, LL,, 9 m.). 
These locations correspond to typical locations of a 
speaker and listener in a large rectangular hall. If we 
consider a scaling factor of about 5, the dimensions of 
this hall would be 115X67X42.2 ft., and our pulse 
would have a corresponding width (at 1/e peak ampli- 
tude) of 8 msec. The scaled carrier frequency would be 
about 700 c.p.s. The resolving time is seen to be shorter 
than the time actually required for the direct pulse to 
reach the receiver, and the character of the response 
shows that a “‘smear” sets in immediately. The actual 
and calculated envelopes are seen to agree well in this 
case, since no degeneracies are present. 

Figures 18 and 19 can be compared qualitatively with 
Fig. 17. These were obtained in a motion picture theater 
of about 250,000 cu. ft. volume. Pulses were produced 
by applying a short pulsed carrier signal to the theater 
loudspeaker. In Fig. 18 the microphone was at a seat 
in the side section of the main floor about two-thirds of 
the way back from the stage, under the balcony. The 
pulse length was 5 msec., and the carrier frequency 
2800 c.p.s. 

In Fig. 19 the microphone was in the center section 
of the main floor, well in front of the balcony and some- 
what off the center line. The pulse length was again 5 
msec. and the carrier frequency was 1500 c.p.s. Qualita- 
tive observations indicated that the “hour glass” bulge 
in Fig. 19 can be correlated subjectively with a “‘slap 
echo” from the hard but irregular rear wall. This bulge 
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Fic. 17. Pulse response for non-degenerate locations of source 
and receiver which were chosen to correspond respectively to the 
positions of a speaker and listener in a large rectangular hall. 
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Fics. 18 AND 19. Typical pulse responses in a motton picture theater of about 250,000 cu. ft. volume. 
Note qualitative similarity with Fig. 17. 


was characteristic of pulse responses at this location for 
nearly all carrier frequencies. A more detailed analysis 
of these and similar photographs is in progress in an 
attempt to correlate various features of the pulse 
response with the results of subjective listening tests. 


CONCLUSIONS 


Most of the theoretical and experimental work pre- 
sented in this paper is directly applicable only to hard- 
walled rectangular rooms. However, as has been 
pointed out, appropriate modifications of the images 
may be made when the walls are not hard and in some 
cases modifications for the first few images do not 
involve too great computational difficulties. 

It appears that fairly detailed experimental and theo- 
retical investigations of the short term transient re- 


sponse (e.g., the first 100-200 msec.) for rooms having 
various shapes and absorptive treatments should prove 
highly instructive. At the same time, the statistical 
methods outlined in this paper can be applied to the 
analysis of the fluctuation characteristics of the long 
term response. 

The experimental results presented show that ap- 
propriate transients can be easily produced and ob- 
served, and, in simple cases at least, correlated with 
boundary conditions. Also it appears that pulse methods 
are readily adaptable to scale model experiments. 
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